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Abstract. We study Langevin dynamics of a driven charged particle in the presence as well as 
in the absence of magnetic field. We discuss the validity of various work fluctuation theorems 
using different model potentials and external drives. We also show that one can generate an orbital 
magnetic moment in a nonequilibrium state which is absent in equilibrium. 
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1. Introduction 

Recent developments in nonequilibrium statistical mechanics has led to the discovery of 
several rigorous theorems for systems far away from equilibrium[l-10]. The fluctua- 
tion theorems describe exact relations for properties (symmetries) of distribution functions 
of various physical quantities such as work, heat, entropy, etc, in the nonequilibrium 
state. The fluctuation relations are statements about the symmetry of the distributions 
around zero and not around maximum of physical quantities . They involve negative 
tails in physical quantities which are usually very rare and are related to transient second 
law violating contributions. These theorems are useful to probe nonequilibrium states 
in nanophysics and biology. In these systems energies involved are typically small and 
hence thermal fluctuations play significant role. In fact, variance in some of the physical 
quantities dominate the mean value rendering these quantities non-self-averaging. Ana- 
lyzing the role of these fluctuations may help in understanding and improving the perfor- 
mance characteristics of engines at nanoscale(e.g, molecular motors). On the application 
side, Jarzynski equality [4] has been used to measure equilibrium free energies (AF)of 
the systems from the statistics of the nonequilibrium work (W) performed. 

There have been an explosion in the number of fluctuation theorems relating various 
physical quantities in the last few years. Some of these theorems have been verified exper- 
imentally on single nanosystems in physical environment where fluctuations play a domi- 
nant role[l 1,12]. In our present work we study some solvable models[8,9,13] which illus- 
trate the Jarzynski equality and related steady state fluctuation theorems[3,8,9,13,14]. We 
have also studied a driven particle in a nonlinear potential numerically to establish steady 
state fluctuation theorem for work . The well- known Bohr-van-Leeuwen theorem 
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states that a classical thermodynamic equilibrium system does not exhibit orbital mag- 
netism[15,16]. However, we show that we can obtain the orbital magnetism (paramag- 
netic/diamagnetic) in driven nonequilibrium systems. 



2. The Model 



We consider the dynamics of a charged(e) Brownian particle in a two dimensional (x — y) 
plane in the presence of a time dependent potential U{=U (x, y, t) ). An external magnetic 
field (B) is along z direction. The particle-environment interactions can be treated via 
Langevin equations [16,17], 

I c I dU 

mx = -^x By- — +£ x (t), (1) 

c ox 



I & I dXJ 

my=- ly +^Bx- — +£„(*), (2) 



where the random force field £ Q (t) is a Gaussian white noise with 

(Za(t)Sf,(i))=D5 a fi5(t-i). (3) 

Here 7 is the friction coefficient and a, /3 = x,y. The consistency conditions for the 
state of equilibrium in the absence of time dependent field relates the prefactor D to 7 as 
D = 2jksT. This problem for time independent potential was considered earlier[15] to 
elucidate the subtle role played by the boundary conditions in the celebrated theorem of 
Bohr-van-Leeuwen on the absence of diamagnetism in classical systems [18]. This, in 
turn, implies that free energy of a system is independent of magnetic field. Correspond- 
ing quantum problem is studied in references [19], with several interesting implications. 
Equations (1) and (2) can be written in the over damped regime as 

\e\B . dU , . . 

7± = -— y-jr + Ut) (4) 

c ox 



t*-^* -£+««> < 5 > 

In our following treatment we consider over damped equations for fluctuation theorems 
and under damped equations for calculating orbital magnetic moment. 

We consider three different protocols for the time dependent potential: (i) Particle in a 
two dimensional harmonic potential the centre of which is dragged with a uniform velocity 
in the diagonal direction in x-y plane. For this case U(x, y, t) = ^k\r — r* | 2 , where 

r is a two dimensional vector (r = xi + yj) and r* (t) = vt(i + j). (ii) U (x, y, t) — 
\k[x 2 + y 2 ) — Ax sin cut, i.e, the particle is subjected to a harmonic ac drive along x 
direction. For case(iii) we consider a one-dimensional problem with nonlinear potential 
U(x) — jax 4 subjected to harmonic drive, i.e, U(x,t) = jax 4 — Ax sin ut. 
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3. Results and Discussions 



Thermodynamic work done on the system for case(i) (or the input energy injected into the 
system ) for an external agent during a time interval t is given by[8,9,20] 

W = -kvf {(x(t')-vt) + (y(t)-vt')}dt, (6) 
Jo 

and for case(ii) and (iii) 

W = -Alu[ cos(ojt)x(t')dt . (7) 
Jo 

It may be emphasized that the thermodynamic work or Jarzynski work is not a mechani- 
cal work[2 1 ] . The thermodynamic work corresponds to the input energy pumped into the 
system by an external time dependent perturbations. It is clear from the above expressions 
that we have to solve the problem for x formally to obtain work distributions. To solve the 
problem analytically for the case(i) we define a new variable z — x + iy (i — \/— T), and 
with the help of the over damped equations (4) and (5), we get 

. -kpz kpg*(t) p£(t) 

z= 1 1 , (8) 

7 7 7 

where p = £(t) - + *&(*). <?*(*) = + and C = ^1. 

The formal solution for equation (8) is given by 

z{t) = z eM--Pt) + - [ dtexp(--p(t-t')){kg*{t)+Z(t')}, (9) 
7 7 Jo 7 

where zq = xo + iyo, and xo and ya are initial co-ordinates of the particle at time 
t = 0. It may be readily noticed from equation (6) that particle co-ordinates at time 
t and consequently work done are linear functionals of Gaussian variables. Hence it 
follows that work distribution[8,9,13] is a Gaussian , which can be completely specified 
by mean(VF) and the variance a 2 = (W 2 ) — (W) 2 . It is straight forward to calculate 
these quantities. For this we refer to [13]. Final result for (W) is given by 

(W) = 2- f v 2 {t - 1(1 - cxp(-fc*t) cos(m)) - ^ sin(m) 

exp(-fc*i)} - "lvHC{\- sin(ftt) exp(-Ft) - ^ 
k k 

(1 - exp(-fc*i) cos(ftt))}, (10) 
where = and k* = ^^y. 

The variance of the work 

a 2 = (W 2 ) - (Wf = HM, (11) 
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Here (3 = jjpr- To obtain the above results (eqns (10) and (1 1)) we have assumed initial 
distribution for the co-ordinates x and y to be equilibrium distribution, P e {x , yo,t) — 
exp[ ~ Wa 2 ° +yg) ] . The full probability distribution P(W) is 

P(W) = 1 P -(w-{w)) 2 /(MW)/0) (U) 

The Jarzynski equality follows immediately namely 

(e-f jw ) = er^ F = 1. (13) 

The above equation implies AF = 0, indicating that the equilibrium free energy 
difference (AF) is independent of magnetic field consistent with Bohr-van-Leeuwen the- 
orem. Jarzynski equality relates nonequilibrium quantities with equilibrium free ener- 
gies. Initially the system is assumed to be in equilibrium defined by a thermodynamic 
parameter A (in our present case centre of the harmonic potential ). The nonequilibrium 
process is obtained by changing the thermodynamic control parameter with a prescribed 
protocol up to time r, where the thermodynamic parameter has value B. The state of the 
system at the end of the protocol is not in equilibrium. This protocol is repeated for large 
number of times. For each realization we get a different quantity W. Using eqn. (1), one 
obtains free energy difference AF = Fb — Fa after evaluating the average (...) over all 
possible realizations. In our present case the free energy is independent of the centre of 
the harmonic oscillator and the applied magnetic field and hence AF = 0. However, it 
may be noticed that the thermodynamic work (eqn. (10)) depends on the magnetic field 

and there is a finite probability of W being negative. The relaxation rate r r (= - ) 
also depends on the magnetic field. In the absence of magnetic field we reproduce the 
results obtained in the references[8,9]. Discussion of the above distribution for P(W) 
in the asymptotic time limit t — > oo it >• r r ) in connection with steady state fluctuation 
theorem and Hatano-Sasa identity is discussed in ref[13]. 

We now turn to case(ii) to examine the steady state fluctuation theorem . In the large 
time regime probability distributions are time periodic with a period (77)- The problem 
being linear we can calculate average work done (W) and variance over a single period ( 
^j) analytically as given by 

2-7T 

(W s ) = lim [(W(t + — )) - (W(t))} (14) 



r^ 7 (fc 2 + W V(l + C 2 )) (15) 



(fc 2 + (1 + C 2 h 2 iu 2 ) 2 - 4k 2 C 2 -f 2 0J 2 



(V s ) = (W 2 ) - (W s ) 2 (16) 

(V.) - |<W a > (17) 
The probability distribution of W s is again Gaussian and satisfies the relation 
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p(-w s ) 



(18) 



The above equation is a statement[8,9,13,14] of steady state fluctuation theo- 
rem[SSFT]. Thus we have shown that work done over a single period in the time asymp- 
totic periodic regime satisfies SSFT We would like to emphasize that the validity of 
SSFT over a single period is restricted only to over damped linear models. In general this 
will not hold true in nonlinear situations. However, we will show later that SSFT holds 
even for nonlinear models if one considers the work done over a large number of periods 
or over a single period , however in the large noise limit. The convergence of SSFT on 
accessible time scales has been discussed in the previous literature [14]. It may also be 
noted that the average work done over a period (en 15) depends on magnetic field. How- 
ever, it is independent of temperature which is again valid for a considered over damped 
linear model only. 

To discuss the validity of SSFT in over damped nonlinear systems we turn to case(iii) 
where the particle in quartic potential is subjected to an ac force . In the absence of 
magnetic field it reduces to a one dimensional problem. Numerical simulations of this 
model was carried out by using Heun's method[22]. To calculate the work done over a 
period (eqn. 7) in the time asymptotic regime we neglect initial transients and work done 
over a period is calculated. To get better statistics we have calculated W s over more than 
1000000 realizations. 

In figs (l-a)-(l-d) we have plotted P(W S ) and P{-W s ) exp(f3W s ) as a function of W s 
over a single period of the ac force for different values of temperature or the noise strength 
(D = 0.02, 0.04, 0.06, 0.6). All the physical parameters are in dimensionless units 
and their values are mentioned in the figure captions. We observe that for small values of 
noise strength SSFT does not hold, i.e., P(W S ) ^ P{— W s )exp{[3W s ). However, only 
in the large noise limit (see fig 1-d) SSFT is indeed satisfied within our numerical accu- 
racy. The large noise (or temperature) limit corresponds to a case where the relaxation 

time(~ vX akgT ^ m tne svstem becomes much less than a given period of ac force. The 
same conclusions can be drawn if one starts in a low noise regime. However, for this 
case one has to evaluate the work done over a large number of cycles such that the relevant 
relaxation time becomes much less than the total time over which the work distribution is 
evaluated. In both the above cases total work can be treated as an addition of independent 
increments (each increment corresponds to work done over a relaxation time). Then the 
central limit theorem leads us to expect that the distribution of work will be Gaussian. And 
this is, indeed, the case. We observe that for D — 0.6, P(W) approaches a Gaussian 
distribution (In the inset of fig (1-d) we have shown a Gaussian fit where the mean and 
variance are calculated from the numerical data). For this distribution the variance V and 
the mean (W) are related by the fluctuation dissipation ratio = 1, so as to satisfy 
SSFT. In our case, for D = 0.6 this ratio is m 0.98. A similar conclusion was arrived 
at recently on work fluctuations in systems exhibiting stochastic resonance[23] . Here it is 
shown that work distribution satisfies SSFT provided one considers work done over a large 
number of periods (low temperature regime mentioned earlier) and in this case distribution 
approaches Gaussian. 

In all the cases studied above we observe that there is always a weight towards negative 
values of W s . The negative values of W s corresponds to the transient second law vio- 
lating hysteresis loops. In the time periodic asymptotic state work done over the periods 
is dissipated into the system as heat[8]. Thus one can identify (W s ) as hysteresis loss 
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(heat) in the medium. However, it may be noted that fluctuations of the work can not be 
identified with heat fluctuations[8]. In fig(2) we have plotted the probability distributions 
of work W, P(W), that of change in internal energy AU, P(AU) and that of heat Q, 
P(Q). All these physical quantities are averaged over a time interval of a single period for 
D = 0.6. From the first law of thermodynamics it follows that W = Q + AU. P(AU) is 
symmetric ({AU) = 0) and the distribution is exponential, i.e., P(AU) ~ e~ ,3 l AC/ l. As 
mentioned earlier the distribution P(W) is Gaussian. At large Q ^> (Q) the distribu- 
tion P(AU) dominates over a Gaussian distribution of P(W) and hence it follows that 
P(Q) ~ e -P\Q-(Q)\ . The presence of this exponential tail in P(Q) , a new extended 
fluctuation theorem for heat is obtained[8]. The consequence of this new theorem for 
heat fluctuations is that the ratio of the probability for the Browinian particle to absorb 
and to supply heat to the environment is much larger than the one corresponding to the 
conventional SSFT for the work. Details of the fluctuation theorem for the heat in our 
nonlinear models will be published elsewhere. We have also observed that the average 
hysteresis loss over a period monotonically decreases with temperature as opposed to the 
case of linear model where it is independent of temperature. 



3 
< 
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W,Q,AU 

Figure 2. Probability distributions for W , Q and AU for temperature D 
other parameters being the same as in fig(l). 



0.6. The 



We finally discuss the problem of orbital magnetism and hysteresis loss in the inertial 
regime. For this we use Langevin equations with inertia (eqns (1) and (2)). The potential 
considered here is U(x, y, t) = \k{x 2 +y 2 ) — Ax sin uit (case-ii). In the time asymptotic 
regime one can readily obtain expressions for the averaged work done over a single period 
{W s } as well as averaged magnetic moment of the system over a period, namely[15], 
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lei lo r t+ « 

(M)= lim[--LL- / (rxv)dt], (19) 
where f and v are the two dimensional position and velocity respectively. We obtain, 

(M) = 2mcVm; ^ V (20) 

where w c = the cyclotron frequency, ft = is the natural frequency of the 

harmonic oscillator and T = The hysteresis loss per period (W s ) is given by 

, s Tr^TcoUn 2 -u; 2 ) 2 + (lu 2 c +T 2 )lo 2 } 

(W s ) = -* ^ ^ ^ J —L (21) 

where 

D = [(cj 2 + T 2 ) 2 ^ + (Q 2 - u; 2 ) 2 {2u; 2 (T 2 - u 2 c ) + (Q 2 - uo 2 ) 2 }} (22) 



<M> 

M 

0-2 - ,' ', O) c =0.5 




Figure 3. The averaged magnetic moment (M) with frequency ui for different values 
of u> c , for k = 1, m = 1, A = 1 and 7 = 1. 



In fig(3) we have plotted dimensionless magnetic moment ( (M) = ( $ )) as a f unc ~ 
tion of dimensionless frequency o> for three different values of cyclotron frequencies. The 
frequencies are scaled with respect to T. In fig(4) we have plotted hysteresis loss in di- 
mensionless form ((W s ) = ( ^^r r )) as a function of frequency uj for different values 
of w c . We have set A, 7 and m to unity. 

From the calculation of (M) we state some noteworthy observations, (i) In a nonequi- 
librium state we obtain magnetic moment, this does not violate Bohr-van Leeuwen theo- 
rem as it is valid only for systems in equilibrium, (ii) (M) goes to zero as u — > 0, this 
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Figure 4. Figure shows the averaged hysteresis loss (W s ) over one period with fre- 
quency lo for different values of lo c . Other parameters are same as in fig(3). 

limit corresponds to equilibrium state, (iii) For small lo, (M) is paramagnetic and 
crosses over to a diamagnetic regime at resonance frequency lo = \J~^ and (M) — » 

as lo — » oo. (iv) In both paramagnetic and diamagnetic regime (M) exhibits a peak 
as a function of w. (v) As B — > 0, (M) — > and again (M) goes to zero in the 
large field limit exhibiting a peak as a function of magnetic field B for a nonzero fixed 
lo. (vi) As a function of friction Y, (M) decreases monotonically and goes to zero 
as T — ► oo. These results are expected on physical grounds. Hysteresis area exhibits 
a double peak behavior as a function of the frequency lo in the inertial regime. This is 
not the case for the over damped motion (see eqn.(15)), in the absence of magnetic field 
(C = 0). These complex structures are attributed to interplay in dynamics with different 
frequencies, namely, cyclotron frequency co c , natural frequency ft along with frequency 
lo of a forcing. 

It may also be noted that in the absence of magnetic field systems do not possess angular 
momentum. However, if we apply two oscillating ac fields along x and y direction 
respectively(in absence of magnetic field) the system can acquire angular momentum as 
the resultant field drags the particle in circular or elliptical orbit depending on the relative 
strength and phase difference between the two perpendicular ac fields. We have also 
solved this problem analytically for (M) as well as the work distributions. Our results 
for work fluctuation in the inertial regime indicate that P{W) is Gaussian and satisfies 
SSFT 



4. Conclusion 

By studying the dynamics of a charged particle in the presence of magnetic field in 
nonequilibrium state we have verified Jarzynski equality. In particular we have shown 
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that nonequilibrium state supports orbital magnetism without violating Bohr-van-Leeuwen 
theorem. Steady state fluctuation theorem for work is discussed for two different model 
systems. In the regime of validity of SSFT we have observed that distribution for work 
approaches a Gaussian distribution. 
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